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It has recently been shown that a Hagedorn phase of string gas cosmology may provide a causal 
mechanism for generating a nearly scale-invariant spectrum of scalar metric fluctuations, without 
the need for an intervening period of de Sitter expansion. A distinctive signature of this structure 
formation scenario would be a slight blue tilt of the spectrum of gravitational waves. In this paper 
we give more details of the computations leading to these results. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

String gas cosmology (SGC) 0,0 (see also for early 
work, jil |^] for recent overviews, and |(| for a critical 
review) is an approach to string cosmology, based on 
adding minimal but crucial inputs from string theory, 
namely new degrees of freedom - string oscillatory and 
winding modes - and new symmetries - T-duality - to the 
now standard hypothesis of a hot and small early universe 
(see for other work on string gas cosmology, and in 
particular f0j for a discussion of the role which D-branes 
play in this scenario). One key aspect of SGC is that the 
temperature cannot exceed a limiting temperature, the 
Hagedorn temperature Tjj Q- This immediately pro- 
vides a qualitative reason which leads us to expect that 
string theory can resolve cosmological singularities 0. 
From the equations of motion of string gas cosmology 
0, 0| it in fact follows that if we follow the radiation- 
dominated Fricdmann-Robertson- Walker (FRW) phase 
of standard cosmology into the past, a smooth transi- 
tion to a quasi-static Hagedorn phase will occur. In this 
phase, the string metric is quasi-static while the dila- 
ton is time-dependent. Reversing the time direction in 
this argument, we can set up the following new cosmo- 
logical scenario Q: the universe starts out in a quasi- 
static Hagedorn phase during which thermal equilibrium 
can be established over a large scale (a scale sufficiently 
large for our current universe to grow out of it following 
the usual non-inflationary cosmological dynamics). The 
quasi-static phase, however, is not a stable fixed point of 
the dynamics, and eventually a smooth transition to a 
radiation-dominated FRW phase will occur, after which 
point the universe evolves as in standard cosmology. This 
transition is connected to the decay of string winding 
modes into string oscillatory degrees of freedom, as de- 
scribed initially in and studied later in more depth 
in [H[ll. 
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Recently 0[T5| it was discovered that, under certain 
assumptions which will be detailed below, string thermo- 
dynamic fluctuations during the Hagedorn phase lead to 
scalar metric fluctuations which are adiabatic and nearly 
scale-invariant at late times, thus providing a simple al- 
ternative to slow-roll inflation for establishing such per- 
turbations. It is to be emphasized that this mechanism 
for the generation of the primordial perturbations is in- 
trinsically stringy - particle thermodynamic fluctuations 
would lead to a spectrum with a large and phenomeno- 
logically unacceptable blue tilt 0] . As discussed in 0] , 
the stringy thermal fluctuations also generate an almost 
scale-invariant spectrum of gravitational waves (tensor 
metric fluctuations). Whereas the spectrum of scalar 
modes has a slight red tilt (like the fluctuations generated 
in most inflationary models), the gravitational wave spec- 
trum has a slight blue tilt. This is a prediction which dis- 
tinguishes our string gas cosmology scenario from stan- 
dard inflationary models. In inflationary models, a slight 
red tilt of the spectrum is expected since the amplitude 
of gravitational waves is set by H, and H is decreasing in 
time, leading to lower values of H when scales with larger 
values of the comoving wavenumber k exit the Hubble ra- 
dius. 

In this paper ; w e review the analyses of our previous 
results [Ij, I15L flq and present more details. The out- 
line of this paper is as follows. We begin with a discus- 
sion of the background string gas cosmology, emphasiz- 
ing the smooth transition between the initial quasi-static 
Hagedorn phase and the later radiation-dominated pe- 
riod. Next, we describe how the scalar and tensor metric 
fluctuations are related to correlation functions of the 
string gas energy-momentum tensor. In Section 4 we 
then discuss the computation of the string gas correla- 
tion functions. Sections 5 and 6 describe the evolution of 
scalar and tensor metric fluctuations, respectively. In the 
final section we discuss our results, point out unresolved 
issues and future directions for research. 

In the following, we assume that our three spatial di- 
mensions are already large (sufficiently large such that 
expansion according to standard cosmology beginning at 
a temperature of the string scale can produce our ob- 
served universe today - if the string scale is of the order 
10 16 GeV, then the initial size needs to be of the order of 
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lmm) during the Hagedorn phase (for a possible mech- 
anism to achieve this see while the extra spatial 
dimensions are confined to the string scale. For a mech- 
anism to achieve the separation into three large and six 
small compact dimensions in the context of SGC see pj 
(see however 0,^^ for caveats), and for a natural dy- 
namical mechanism arising from SGC to stabilize all of 
the moduli associated with the extra spatial dimensions 

see ei in m m mm (see also 241). 

To be specific, we take our three dimensions to be 
toroidal. The existence of one cycles results in the stabil- 
ity of string winding modes - and this is a key ingredient 
in our calculations. 



II. THE BACKGROUND 

String gas cosmology 0, is based on coupling an 
ideal gas of strings to a dilaton-gravity background, and 
on focusing on the consequences which follow from the ex- 
istence of new stringy degrees of freedom and new stringy 
symmetries, namely degrees of freedom and symmetries 
which are not present in point particle theories. 

In the case of a homogeneous and isotropic background 
given by 

ds 2 = dt 2 - a(t) 2 dx 2 , (1) 

the three resulting equations of motion of dilaton-gravity 
(the generalization of the two Friedmann equations plus 
the equation for the dilaton) in the string frame are Q 
(see also \W^ ) 

-{d)\ 2 +cp 2 = e v E (2) 
A-^A = l -e*P (3) 

(p-{d)\ 2 = ±e*E, (4) 

where E and P denote the total energy and pressure, 
respectively, d is the number of spatial dimensions, and 
we have introduced the logarithm of the scale factor 

A(i) = In [a(i)], (5) 

and the rescaled dilaton 

if = 2<j)-{d)\. (6) 

If we take the spatial topology to be that of a torus 
with radius R, the new stringy degrees of freedom are 
string winding modes whose energies are quantized in 
units of R, and oscillatory modes, whose energies are 
independent of R. The number of oscillatory modes of a 
string increases exponentially with energy, thus leading 
to a maximal temperature (the Hagedorn temperature 
Th 0) for a gas of strings. 

Since the energies of the momentum modes of a string 
are quantized in units of l 2 /R, where l s is the string 



length, the spectrum of mass states of the free string 
is symmetric under the interchange R — > l 2 /R, the T- 
duality transformation. String vertex operators are also 
invariant under this transformation, and assuming that 
the symmetry extends to non-perturbative string theory 
leads to the existence of D-branes as new fundamental 
objects in string theory [25| . 

In thermal equilibrium at the string scale (R ~ l s ), 
the self-dual radius, the number of winding and momen- 
tum modes will be equal. Since winding and momentum 
modes give an opposite contribution to the pressure, the 
pressure of the string gas in thermal equilibrium at the 
self-dual radius will vanish. From the above equations of 
motion it then follows that a static phase A = will be a 
fixed point of the dynamical system. This phase is called 
the "Hagedorn phase" . 

On the other hand, for large values of R in thermal 
equilibrium the energy will be exclusively in momentum 
modes. These act as usual radiation. Inserting the radia- 
tive equation of state into the above equations 10 -0J it 
follows that the source in the dilaton equation of motion 

4> + (d)\<P - 20 2 = -\e*E[l - (d)w] , (7) 

where w = P/E is the equation of state parameter, van- 
ishes and the dilaton approaches a constant as a con- 
sequence of Hubble damping. Consequently, the scale 
factor expands as in the usual radiation-dominated uni- 
verse. 

The transition between the Hagedorn phase and the 
radiation-dominated phase with fixed dilaton is achieved 
via the annihilation of winding modes, as studied in detail 
in The main point is that, starting in a Hagedorn 
phase, there will be a smooth transition to the radiation- 
dominated phase of standard cosmology with fixed dila- 
ton. It is in this cosmological background that we will 
study the generation of fluctuations. 

It is instructive to compare the background evolution 
of string gas cosmology with the background of inflation- 
ary cosmology |26l l27j |. the current paradigm of early 
Universe evolution. Figure 1 is a sketch of the space- 
time evolution in string gas cosmology. For times t < tu, 
we are in the quasi-static Hagedorn phase, for t > tn , we 
have the radiation-dominated period of standard cosmol- 
ogy. To understand why string gas cosmology can lead to 
a causal mechanism of structure formation, we must com- 
pare the evolution of the physical wavelength correspond- 
ing to a fixed comoving scale (fluctuations in early uni- 
verse cosmology correspond to waves with a fixed wave- 
length in comoving coordinates) with that of the Hubble 
radius H^ 1 ^), where H(t) is the expansion rate. The 
Hubble radius separates scales on which fluctuations os- 
cillate (wavelengths smaller than the Hubble radius) from 
wavelengths where the fluctuations are frozen in and can- 
not be effected by microphysics. Causal microphysical 
processes can gene rate fluctuations only on sub-Hubble 
scales (see e.g. [28|| for a concise overview of the theory of 



FIG. 1: Space-time diagram (sketch) showing the evolution 
of fixed comoving scales in string gas cosmology. The verti- 
cal axis is time, the horizontal axis is physical distance. The 
Hagedorn phase ends at the time tR and is followed by the 
radiation-dominated phase of standard cosmology. The solid 
curve represents the Hubble radius H~ x which is cosmological 
during the quasi-static Hagedorn phase, shrinks abruptly to a 
microphysical scale at tR and then increases linearly in time 
for t > tR. Fixed comoving scales (the dotted lines labeled by 
ki and k%) which are currently probed in cosmological obser- 
vations have wavelengths which are smaller than the Hubble 
radius during the Hagedorn phase. They exit the Hubble ra- 
dius at times ti(k) just prior to tR, and propagate with a 
wavelength larger than the Hubble radius until they reenter 
the Hubble radius at times tf(k). 

cosmological perturbations and 29] for a comprehensive 
review) . 

In string gas cosmology, the Hubble radius is infinite 
deep in the Hagedorn phase. As the universe starts ex- 
panding near tR, the Hubble radius rapidly decreases to 
a microscopic value (set by the temperature at t = tR 
which will be close to the Hagedorn temperature), be- 
fore turning around and increasing linearly in the post- 
Hagedorn phase. The physical wavelength corresponding 
to a fixed comoving scale, on the other hand, is con- 
stant during the Hagedorn era. Thus, all scales on which 
current experiments measure fluctuations are sub-Hubble 
deep in the Hagedorn phase. In the radiation period, the 
physical wavelength of a perturbation mode grows only 
as \/t. Thus, at a late time tf(k) the fluctuation mode 
will re-enter the Hubble radius, leading to the perturba- 
tions which are observed today. 



FIG. 2: Corresponding space-time diagram (sketch) for infla- 
tionary cosmology. The period of exponential expansion is for 
t < tR. At the time of reheating t — tR the universe makes 
a transition to a radiation-dominated phase. As in Figure 1, 
the vertical axis is time and the horizontal axis is physical dis- 
tance. The solid curve represents the Hubble radius. Fixed 
comoving scales (the dotted lines labeled by fci and fe) which 
are currently probed in cosmological observations have wave- 
lengths which start out smaller than the Hubble radius at the 
beginning of the inflationary phase. They exit the Hubble 
radius at times U(k) and propagate with a wavelength larger 
than the Hubble radius until they reenter the Hubble radius 
at times tf(k). 



In contrast, in inflationary cosmology (Figure 2) the 
Hubble radius is constant during inflation (t < tji, where 
here tR is the time of inflationary reheating) , whereas the 
physical wavelength corresponding to a fixed comoving 
scale expands exponentially. Thus, as long as the period 
of inflation is sufficiently long, all scales of interest for 
current cosmological observations are sub-Hubble at the 
beginning of inflation. 

In spite of the fact that both in inflationary cosmology 
and in string gas cosmology, scales are sub-Hubble dur- 
ing the early stages, the actual generation mechanism 
for fluctuations is completely different. In inflationary 
cosmology, any thermal fluctuations present before the 
onset of inflation are red-shifted away, leaving us with a 
quantum vacuum state, whereas in the quasi-static Hage- 
dorn phase of string gas cosmology matter is in a ther- 
mal state. Hence, whereas in inflationary cosmology the 
fluctuations originate as quantum vacuum perturbations 
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30, 31], in string gas cosmology the inhomogencities are 
created by the thermal fluctuations of the string gas. 

As we have shown in [Til [TH Il6[ . string thermody- 
namics in the Hagedorn phase of string gas cosmology 
yields an almost scale-invariant spectrum of both scalar 
and tensor modes. This is an intrinsically stringy result: 
thermal fluctuations of a gas of particles would lead to a 
very different (namely very "blue" spectrum). Since the 
primordial perturbations in our scenario are of thermal 
origin (and there are no non-vanishing chemical poten- 
tials), they will be adiabatic. 

In the following, we will discuss the derivation of these 
results in more details than discussed in our initial papers 

MM 



III. SCALAR AND TENSOR METRIC 
FLUCTUATIONS 



In this section, we show how the scalar and tensor met- 
ric fluctuations can be extracted from knowledge of the 
energy-momentum tensor of the string gas. 

Working in conformal time rj (defined via dt = a(t)drj), 
the metric of a homogeneous and isotropic background 
space-time perturbed by linear scalar metric fluctuations 
and gravitational waves can be written in the form 



ds 2 = a 2 (ry) {(1 + ^dry 2 - [(1 - + /i^cfoW} . 

(8) 

Here, <f> and \t describe the scalar metric fluctuations. 
Both are functions of space and time. The tensor hij is 
transverse and traceless and contains the two polarization 
states of the gravitational waves. In the above, we have 
fixed the coordinate freedom by working in the so-called 
"longitudinal" gauge in which the scalar metric fluctu- 
ation is diagonal. Note that to linear order in the am- 
plitude of the fluctuations, scalar and tensor modes de- 
couple, and the tensor modes are gauge- invariant. Note 
also that in the absence of non-adiabatic stress, the two 
degrees of freedom $ and \I/ for scalar metric fluctuations 
coincide (see |29j for a comprehensive review of the theory 
of cosmological perturbations and [28| for a pedagogical 
introduction). 

Inserting the above metric into the Einstein equations, 
subtracting the background terms and truncating the 
perturbative expansion at linear order leads to the fol- 
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where D = $ — \P, Tl = a' /a, a prime denotes the deriva- 
tive with respect to conformal time rj, and G is Newton's 
gravitational constant. 

In the Hagedorn phase, these equations simplify sub- 
stantially and allow us to extract the scalar and tensor 
metric fluctuations individually. First of all, since there 
is no anisotropic stress at late times, D — and hence 
<J> = . Replacing comoving by physical coordinates, we 
obtain from the 00 equation 



V 2 $ = 4ttG5T° 



and from the i ^ j equation 



V 2 ^ 



-AttGST 1 



(10) 



(11) 



The above equations (|10fl and i|ll|) allow us to compute 
the power spectrum of scalar and tensor metric fluctu- 
ations in terms of correlation functions of the energy- 
momentum tensor. Since the metric perturbations are 
small in amplitude we can consistently work in Fourier 
space. 

Our approximation scheme for computing the cosmo- 
logical perturbations and gravitational wave spectra from 
string gas cosmology is as follows (the analysis is simi- 
lar to how the calculations were performed in [321 33] in 
the case of inflationary cosmology) . For a fixed comoving 
scale k we follow the matter fluctuations until the time 
U(k) shortly before the end of the Hagedorn phase when 
the scale exits the Hubble radius. At that time, we use 
(|10p and Ijllfl to compute the values of $(fc) and h(k) (h 
is the amplitude of the gravitational wave tensor), and, 
since we are working in Fourier space, we replace the V 2 
operator by k 2 . 

The expectation value of the metric fluctuations is 
thus given in terms of correlation functions of the string 
energy-momentum tensor. Specifically, 

(|$(fc)| 2 > = 16^ 2 G 2 fc^ 4 ( ( 5T (fc) ( 5T o(fc)) , (12) 

where the pointed brackets indicate expectation values, 
and 



(\h(k)\ 2 ) = 16TT 2 G 2 k^(5T i J (k)ST i J (k)) 



(13) 
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where on the right hand side of 113fl we mean the average 
over the correlation functions with i ^= j. 

After the modes exit the Hubble radius, we follow the 
evolution of both scalar and tensor fluctuations by means 
of the usual equations of cosmological perturbation the- 
ory. Apart from small corrections which arise in the time 
interval t G \ti(k), to], this technique is justified since the 
dilaton is fixed after the end of the Hagedorn phase. 

In the following section, we describe in a bit more detail 
than in our previous letters how to compute the string 
matter correlation functions from string thermodynam- 
ics. After that, we return to the computation of the 
spectrum of scalar and tensor modes. 



and hence 



(T"„) = 2- 



t x S In Z 



^ fig 



(20) 



To extract the fluctuation tensor of T M „ for long wave- 
length modes, we take one additional variational deriva- 
tive of (J20J, using [(TSj) to obtain 



g aS d\nZ^ 




IV. ENERGY-MOMENTUM TENSOR 
CORRELATION FUNCTIONS FOR CLOSED 
STRINGS 

In the following we will denote the space-time metric 
by g^. Unlike what was done in |16| . we will in this pa- 
per be working with a Minkowski signature metric. The 
starting point for our considerations is the free energy F 
of a string gas in thermal equilibrium 



As we saw in the previous section, the scalar metric 
fluctuations are determined by the energy density corre- 
lation function 



cVo = (fip 2 ) = (?) (pY 



(22) 



We will read off the result from the expression (12111 evalu- 
ated for fi = h , = a = X = 0. The derivative with respect 
to goo can be expressed in terms of the derivative with 
respect to (3. After a couple of steps of algebra we obtain 



-~lnZ, 



(14) 



where (3 is the inverse temperature and the canonical 
partition function Z is given by 



Z = 



(15) 



where the sum runs over the states s of the string gas, 
and H(s) is the energy of the state. The action S of the 
string gas is given in terms of the free energy F via 



S = J dty/-gooF[gij,/3} 



(16) 



Note that the free energy depends on the spatial compo- 
nents of the metric because the energy of the individual 
string states does. The energy-momentum tensor T^ v of 
the string gas is determined by varying the action with 
respect to the metric: 
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(17) 



Consider now the thermal expectation value 



(r"„) = ^r"„( s K^ ff ' s ', (is) 



where T M „(s) and H(s) are the energy momentum tensor 
and the energy of the state labeled by s, respectively. 
Making use of (|17[) and l|16|l we immediately find that 



T\{a) = 2- 
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where 



C v = (dE/dT)\ v 



is the specific heat, and 



-MS) 



(23) 



(24) 



(25) 



is the internal energy. Also, V = R 3 is the volume of the 
three compact but large spatial dimensions. 

The gravitational waves are determined by the off- 
diagonal spatial components of the correlation function 
tensor, i.e. 



C 



j j 



= (ST'/) = (TV) - <r,) : 



(26) 



with i ^= j. 

Our aim is to calculate the fluctuations of the energy- 
momentum tensor on various length scales R. For each 
value of R, we will consider string thermodynamics in 
a box in which all edge lengths are R. From 121|) it is 
obvious that in order to have non-vanishing off-diagonal 
spatial correlation functions, we must consisted a twisted 
torus. Let us focus on the x — y component of the corre- 
lation function. We will consider the spatial part of the 
metric to be 



ds 2 



R 2 del + 2eR z d9 x dO v + R z d& 2 



(27) 



with e <C 1. The spatial coordinates Qi run over a fixed 
interval, e.g. [0, 2tv]), The generalization of the spatial 
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part of the metric to three dimensions is obvious. At the 
end of the computations, we will set e = 0. 

From the form of l|21l) . it follows that all space-space 
correlation function tensor elements are of the same order 
of the magnitude, namely 



3 3 



1 



d 



1 dF 



l3R 3 d\nR \ i? 3 <91ni? 
where the string pressure is given by 

--(—) =r(- 
V\d\nR \dV 



1 dp 



P 



(28) 



(29) 



In the following, we will compute the two correlation 
functions (|23|) and (|28J) using tools from string statistical 



mechanics. Specifically, we will be following the discus- 
sion in (see also [Hi EE El H| 
is the formula 



). The starting point 



S(E,R) = \nn(E,R) 



(30) 



for the entropy in terms of £l(E, R), the density of states. 
The density of states of a gas of closed strings on a large 
three-dimensional torus (with the radii of all internal di- 
mensions at the string scale) was calculated in Q (see 
also 0) and is given by 



n(E,R) ~ (3 H e PH 



E+n H V 



[l + Sn (1) (E,R)} 



(31) 



where <5f2(i) comes from the contribution to the density 
of states (when writing the density of states as a Laplace 
transform of Z((3), which involves integration over (3) 
from the closest singularity point Pi to Ph = (1/Tff) 
in the com plex (3 plane. Note that (3\ < Ph, and P\ is 
real. From (IHk^H we have 



SSl w (E,R) 



_ (PhE) 5 -(Ph-p^e-phV) 
5! 



(32) 



In the above, nn is a (constant) number density of order 
l~ 3 and ph is the 'Hagedorn Energy density' of the order 
l~ 4 , and 

Pu ~ 01 \ (R 2 /l s ) , for R^l s . (33) 
To ensure the validity of Eq. 1)3 1|1 we demand that 



sn w < 1 



(34) 



by assuming p = (E/V) ^> ph- 

Combining the above results, we find that the entropy 
of the string gas in the Hagedorn phase is given by 



S(E, R) ~ (3 H E + n H V + In [l + SQ {1) ] 



(35) 



and therefore the temperature T(E, R) = \{dS/dE)v] 1 
will be 



T 



T, 



H 



dSn {1) /dE 

i + sn {1) 

0H-01 



Ph 



(36) 



In the above, we have dropped a term which is negligible 
since E(/3h — Pi) ^ 1 (see (GU). Using this relation, we 
can express tffim in terms of T and R via 



i 3 s sn (1 



(i) 

In addition, we find 

E ~ C 3 i? 2 ln 



I? 



T 
Th~ 



i? 2 (l - T/T H ) 



(37) 



(38) 



Note that to ensure that |<5^(i)| -C 1 and E ^> pnR 3 ! 
one should demand 



(1 - T/T H )R 2 l- 2 < 1 . 



(39) 



The results (fTTH|> and (|3*7jl now allow us to compute the 
correlation functions (|23|l and 1)28(1 . We first compute the 
energy correlation function l|23|) . Making use of (|38H . it 
follows from that 



R 2 ni 



T(1-T/T H ) 



from which we get 



C 00 







(S P 2 



T 



ll{l-T/T H )R 4 



(40) 



(41) 



Note that the factor (1 — T/Tjj) in the denominator is 
responsible for giving the spectrum a slight red tilt. It 
comes from the differentiation with respect to T. 

Next we evaluate ()29|l . From the definition of the pres- 
sure it follows that (to linear order in 50(1)) 



d 

p = n H T + T—6fl {1) 



(42) 



where the final partial derivative is to be taken at con- 
stant energy. In taking this partial derivative, we insert 
the expression Ipf7|l for and must keep careful ac- 

count of the fact that Ph — Pi depends on the radius R. 
In evaluating the resulting terms, we keep only the one 
which dominates at high energy density. It is the term 
which comes from differentiating the factor Pjj — Pi . This 
differentiation brings down a factor of E, which is then 
substituted by means of (f 381) , thus introducing a logarith- 
mic factor in the final result for the pressure. We obtain 



P{E, R) ~ n H T H - -jjj^ In 



R?(1-T/T H ) 



which immediately yields 

O 4 4 



T(l - T/T H ) ^ 



11 R 4 



;i - T/T H ) 



(43) 



(44) 



Note that since no temperature derivative is taken, the 
factor (1 — T/Tjj) remains in the numerator. This is 
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one of the two facts which leads to the slight blue tilt of 
the spectrum of gravitational waves. The second factor 
contributing to the slight blue tilt is the explicit factor of 
R 2 in the logarithm. Because of l|39[) . we are on the large 
k side of the zero of the logarithm. Hence, the greater the 
value of k, the larger the absolute value of the logarithmic 
factor. 



V. POWER SPECTRUM OF SCALAR METRIC 
FLUCTUATIONS 

Let us recall the philosophy behind our computation 
of the spectrum of scalar metric fluctuations. For a fixed 
comoving scale labeled by its wavenumber k, we follow 
the energy-momentum tensor of the string matter fluctu- 
ations until the time U(k) when the scale exits the Hub- 
ble radius. At that point, we use the Einstein constraint 
equation 112fl to determine the corresponding metric fluc- 
tuations. On super- Hubble scales for t > U(k) we use the 
usual equations for cosmological perturbations to evolve 
the fluctuations. 

The dimensionless power spectrum of the cosmological 
perturbations (scalar metric fluctuations) is defined by 
(see e.g. [33) 



P 9 (k) = fc J |$(fc)| 



(45) 



where we are using the convention of defining the Fourier 
transform of a function f(x) including a factor of the 
square root of the volume of space. Thus, the dimension 
of $(fc) is k~ 3 / 2 . 

Making use of (|12[1 to replace the expectation value of 
|$(fc)| 2 in terms of the correlation function of the energy 
density, we obtain 



P»(k) = l6Tr 2 G 2 k- x (\5p(k)f 



(46) 



The Fourier space expectation value of Sp is related to 
the position space root mean square density contrast SM 
in a region of radius R — k^ 1 via 



((6M) 2 ) R = k- 3 (\S P (k)\ 2 ), 

and hence (|46|l becomes 

P$(fc) = l6ir 2 G 2 k 2 ((6M) 2 ) R 
= 16tt 2 G 2 k- i {(Sp) 2 ) B 



(47) 



(48) 



where in the second step we have replaced the mass fluc- 
tuation by the position-space density fluctuation. Ac- 
cording to (|23[) . the density correlation function is given 
by the specific heat via T 2 R~ 6 Cy, and hence (|4*5|> be- 
comes 

P 9 (k) = l6ir 2 G 2 k 2 T 2 C v ■ (49) 

Inserting the expression from (|40() for the specific heat of 
a string gas on a scale R yields our final result 



P$(fc) = 16tt 2 G 



T 



1 



II 1 - T/T H 



(50) 



for the power spectrum of cosmological fluctuations. 

During the Hagedorn phase, the temperature T is close 
to the Hagedorn temperature Tjj , which in turn is given 
by the string scale. Thus, to a first approximation, when 
evaluating the spectrum l|50l) at the time ti(k), the tem- 
perature T(ti(k)) is independent of k. Thus, the spec- 
trum of scalar metric fluctuations is scale-invariant, with 
an amplitude given by 



1 



1-T(ti(k))/T H 



(51) 



The last factor yields a small red tilt of the spectrum 
since T(ti(k)) is a slowly decreasing function of k (modes 
with larger values of k are exiting the Hubble radius later, 
i.e. at a slightly lower value of T(ti(k))). 

As long as the time-dependence of the dilaton is negli- 
gible, then on super-Hubble scales the quantity £ which 
corresponds to the curvature perturbation in comoving 
gauge is constant (modulo effects which come from the 
decaying mode of <I> at U(k) 0]). The quantity ( is 
defined by |H 



c = $ + 



2 ( J ff' 1 $ + $) 



1 



(52) 



where w — p/p is the equation of state parameter. 

Unlike in inflationary cosmology, where the factor 1+w 
changes by many orders of magnitude during reheating, 
in string gas cosmology 1+w makes a very mild transition 
between an initial value of (1 + u>)|j = 1 and a final value 
of (1 + w)\f = 4/3 between the Hagedorn phase and the 
radiation-dominated phase. Thus, up to a factor of order 
unity, the conservation of £ on super-Hubble scales in our 
case leads to the conservation of $. 

It is important to realize, however, that in spite of the 
fact that the amplitude of $ is constant on super-Hubble 
scales, the fluctuation mode is evolving non-trivially - it 
is being highly squeezed by the expansion of the back- 
ground. One way to see this is to realize that when the 
fluctuations cross the Hubble radius, $ is comparable to 
k<S> (since the fluctuations are thermal in origin). How- 
ever, for t > ti(k) the amplitude of <E> decreases pro- 
portionally to a _1 (£). The fluctuation mode approaches 
the form of a standing mode at rest. This squeezing 
is analogous to what happens in inflationary cosmology. 
The only difference is that in inflationary cosmology, the 
squeezing takes place both during and after inflation, 
whereas in string gas cosmology the squeezing happens 
only during the period of radiation-domination. 

Another way to understand the squeezing of cosmo- 
logical fluctuations is in terms of canonical variables 
(variables characterizing the cosmological fluctuations 
in terms of which the action for the perturbations has 
canonical kinetic term). One of the key results of the 
quantum theory of cosmological fluctuations [3 |2^ is 
that the variable v in terms of which the action for fluc- 
tuations (for ideal gas matter) has the canonical form 
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S = 



d A x 



C s V,iV,i 



is the so-called Sasaki-Mukhanov variable I42t 14 



(tp v - 2z$) 



(53) 



(54) 



where <p v is the velocity potential of the ideal fluid, c s is 
the speed of sound, and z is a function of the background 
given by 



z(r)) = a{rj) 
oc a(ji) . 



yn 2 - w 

c s H 



(55) US: 



The action i|53|l is that of a free scalar field with a time- 
dependent negative square mass. The mass term domi- 
nates on super-Hubble scales and leads to the squeezing 
of the wave function for v, squeezing which is propor- 
tional to a. 

The squeezing of fluctuations is crucial in order to 
obtain the "acoustic" oscillations in the angular power 
spectrum of the cosmic microwave background (CMB) 
anisotropies 0] . If all modes enter the Hubble radius at 
late times as standing waves and then begin to oscillate, 
then at the time of last scattering t^g the wave which is 
entering the Hubble radius at that time will lead to max- 
imal redshift fluctuations, whereas a mode which has en- 
tered the Hubble radius earlier and has had time to per- 
form 1/4 oscillation will have zero amplitude and max- 
imal velocity at tj,s, and on the corresponding angular 
scale the redshift contribution to the CMB anisotropies 
will vanish. Models without squeezing of fluctuations on 
super-Hubble scales do not yield this pattern of acoustic 
oscillations, even if they, such as in the case of topologi- 
cal defect models of structure formation , generate a 
scale- invariant spectrum of fluctuations. 



VI. POWER SPECTRUM OF TENSOR METRIC 
FLUCTUATIONS 

In this section we estimate the dimensionless power 
spectrum of gravitational waves. First, we make some 
general comments. In slow-roll inflation, to leading order 
in perturbation theory matter fluctuations do not couple 
to tensor modes. This is due to the fact that the mat- 
ter background field is slowly evolving in time and the 
leading order gravitational fluctuations are linear in the 
matter fluctuations. In our case, the background is not 
evolving (at least at the level of our computations), and 
hence the dominant metric fluctuations are quadratic in 
the matter field fluctuations. At this level, matter fluctu- 
ations induce both scalar and tensor metric fluctuations. 
Based on this consideration we expect that in our string 



gas cosmology scenario, the ratio of tensor to scalar met- 
ric fluctuations will be larger than in simple slow-roll in- 
flationary models. 

The method for calculating the spectrum of gravita- 
tional waves is similar to the procedure outlined in the 
last section for scalar metric fluctuations. For a mode 
with fixed comoving wavenumber k, we compute the cor- 
relation function of the off-diagonal spatial elements of 
the string gas energy- momentum tensor at the time U(k) 
when the mode exits the Hubble radius and use (|T3|) to 
infer the amplitude of the power spectrum of gravita- 
tional waves at that time. We then follow the evolution of 
the gravitational wave power spectrum on super-Hubble 
scales for t > U(k) using the equations of general rela- 
tivistic perturbation theory. 

The power spectrum of the tensor modes is given by 



P h (k) = 16Tr 2 G 2 k- 4 k i (6T l ] (k)5T l j (k)} (56) 

for i 7^ j. Note that the A; 3 factor multiplying the mo- 
mentum space correlation function of T % j gives the po- 
sition space correlation function, namely the root mean 
square fluctuation of T % j in a region of radius R = fc _1 
(the reader who is skeptical about this point is invited to 
check that the dimensions work out properly). Thus, 

167r 2 G 2 fc- 4 C l 



Ph{k) 



(57) 



The correlation function C l j % j on the right hand side of 
the above equation was computed in Section 4. Inserting 
the result (|44|l into (|56|l we obtain 



P h {k) ~ 16ir 2 G 2 -r{l-T/T H )\n 



^- T I T h) 



which, for temperatures close to the Hagedorn value re- 
duces to 



Ph(k) 



l_Pl 
I, 



(l-T/T H )ln 2 



1 



l 2 k 2 



(1 - T/T H ) 



(59) 

This shows that the spectrum of tensor modes is - to 
a first approximation, namely neglecting the logarith- 
mic factor and neglecting the fc-dependence of T(ti(k)) 
- scale- invariant. The corrections to scale-invariance will 
be discussed at the end of this section. 

On super-Hubble scales, the amplitude h of the grav- 
itational waves is constant. The wave oscillations freeze 
out when the wavelength of the mode crosses the Hub- 
ble radius. As in the case of scalar metric fluctuations, 
the waves are squeezed. Whereas the wave amplitude re- 
mains constant, its time derivative decreases. Another 
way to see this squeezing is to change variables to one 



^(ry,x) = a(i])h(r],x.) 



(60) 



in terms of which the action has canonical kinetic term. 
The action in terms of ij) becomes 



S 



l - J d^x (V - titi + ^V) (61) 
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from which it immediately follows that on super-Hubble 
scales ip ~ a. This is the squeezing for gravitational 
waves. 

Since there is no fc-dependence in the squeezing fac- 
tor, the scale-invariance of the spectrum at the end of 
the Hagedorn phase will lead to a scale-invariance of the 
spectrum at late times. 

Note that in the case of string gas cosmology, the 
squeezing factor 2(77) does not differ substantially from 
the squeezing factor a{rf) for gravitational waves. In 
the case of inflationary cosmology, z{rf) and a(rj) dif- 
fer greatly during reheating, leading to a much larger 
squeezing for scalar metric fluctuations, and hence to a 
suppressed tensor to scalar ratio of fluctuations. In the 
case of string gas cosmology, there is no difference in 
squeezing between the scalar and the tensor modes. This 
supports our expectation that the tensor to scalar ratio 
will be larger in our scenario than in typical single field 
inflationary models. We do find a relative suppression of 
tensor modes, but its origin comes from string thermo- 
dynamics in the Hagedorn phase. 

Let us return to the discussion of the spectrum of 
gravitational waves. The result for the power spectrum 
is given in (|59|l . and we mentioned that to a first ap- 
proximation this corresponds to a scale-invariant spec- 
trum. As realized in |16| . the logarithmic term and the 
/c-dependence of T(ti(k)) both lead to a small blue-tilt of 
the spectrum. This feature is characteristic of our sce- 
nario and cannot be reproduced in slow-roll inflationary 
models. 

To study the tilt of the tensor spectrum, we first have 
to keep in mind that our calculations are only valid in 
the range (|39[) . i.e. to the large k side of the zero of the 
logarithm. Thus, in the range of validity of our analysis, 
the logarithmic factor contributes an explicit blue tilt of 
the spectrum. The second source of a blue tilt is the fac- 
tor 1 — T(ti(k)) /Tjj multiplying the logarithmic term in 
(|59|l . Since modes with larger values of k exit the Hub- 
ble radius at slightly later times U(k), when the temper- 
ature T(ti(k)) is slightly lower, the factor will be larger. 
Figure 3 shows the spectrum in the example of an in- 
stantaneous transition between the Hagedorn phase and 
the radiation-dominated period (i.e. no fc-dependence of 
T( tl (k))). 

Note that in the result for the amplitude of the scalar 
metric fluctuations, the prefactor 1 — T(ti{k)) /Th ap- 
pears in the denominator and hence leads to a slight red 
tilt of the spectrum. In addition, the logarithmic factor 
is absent. 

A heuristic way of understanding the origin of the 
slight blue tilt in the spectrum of tensor modes is as fol- 
lows. The closer we get to the Hagedorn temperature, 
the more the thermal bath is dominated by long string 
states, and thus the smaller the pressure will be com- 
pared to the pressure of a pure radiation bath. Since the 
pressure terms (strictly speaking the anisotropic pressure 
terms) in the energy-momentum tensor are responsible 
for the tensor modes, we conclude that the smaller the 




FIG. 3: Power spectrum of gravitational waves in string gas 
cosmology in the case of an instantaneous transition between 
the Hagedorn phase and the radiation-dominated period. In 
this approximation, increasing the factor 1 — T /Th is equiva- 
lent to creasing R. The range of values of R where our analysis 
applies to the left of the zero of the curve at the value (l s /R) 2 . 
The increase of the power spectrum (the vertical axis) going 
towards the left shows that the spectrum has a blue tilt. 



value of the wavenumber k (and thus the higher the tem- 
perature T(ti(k)) when the mode exits the Hubble radius, 
the lower the amplitude of the tensor modes. In contrast, 
the scalar modes are determined by the energy density, 
which increases at T{U(k)) as k decreases, leading to a 
slight red tilt. 



VII. DISCUSSION 

In this paper we have presented more details of the 
derivation of the power spectra for scalar and tensor met- 
ric fluctuations. Here we summarize the main results. 

First of all, we have found that our string gas cos- 
mology background leads to a new way of obtaining a 
roughly scale-invariant spectrum of cosmological pertur- 
bations. Our mechanism relies on string theory in two 
crucial ways. Firstly, the cosmological background is 
the result of a new symmetry of string theory which is 
not present in ordinary quantum field theory, namely T- 
duality. The background equations of motion which are 
consistent with the T-duality symmetry, namely those 
of dilaton gravity, predict that the radiation-dominated 
phase of standard cosmology was preceded by a phase 
- the Hagedorn phase - during which the scale factor 
in the string frame is static, the pressure vanishes, and 
the temperature is close to the limiting temperature, the 
Hagedorn temperature. As we have shown, string ther- 
modynamical fluctuations in the Hagedorn phase lead to 
an almost scale-invariant spectrum of cosmological per- 
turbations. This is the second point where string theory 
plays a crucial role. Point particle fluctuations would 
have given a completely different spectrum. 

We have shown that the spectrum of the scalar metric 
fluctuations has a slight red tilt, whose value depends on 
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how fast the transition between the Hagedorn phase and 
the radiation-dominated period occurs. 

Our model makes a unique prediction: the spectrum 
of gravitational waves has a slight blue tilt, unlike the 
slight red tilt which is predicted in scalar field-driven in- 
flationary models. It would be interesting to investigate 
how hard it would be to observationally identify this sig- 
nature. One factor which will help in this regard is the 
fact that the tensor to scalar ratio r in our scenario is 
not suppressed by slow-roll parameters. As can be seen 
immediately by comparing our results for the amplitude 
of the tensor and scalar power spectra, respectively, the 
ratio r at a scale k is given by 



r ~ {l~T(U{k))/T H f\n 2 



1 



l 2 k 2 



(l-T(U(k))/T H ) 



In principle (if the dynamical evolution from the Hage- 
dorn phase to the radiation-dominated FRW phase were 
under complete analytical control) this quantity would 
be calculable. If the string length were known, the fac- 
tor (1 — T/Th) could be determined from the normaliza- 
tion of the power spectrum of scalar metric fluctuations. 
Since the string length is expected to be a couple of or- 
ders larger than the Planck length, the above factor does 
not need to be very small. Thus, generically we seem to 
predict a ratio r larger than in simple slow-roll inflation- 
ary models. 

We close by mentioning some limitations of our anal- 
ysis. Since the dilaton is evolving until the beginning 
of the radiation phase, the gravitational perturbations 



should be determined and evolved using the complete 
dilaton gravity background, not using the equations of 
general relativity. The time dependence of the dilaton 
will thus lead to correction terms. We are currently in 
the process of calculating the magnitude of these effects. 
This calculation will also be crucial in order to give a 
correct estimate of the amplitude of the slight tilt in our 
spectra. These calculations are also currently in progress. 

We should also emphasize that, at least in the frame- 
work presented in this paper, our cosmological scenario 
does not solve the entropy and flatness problems of stan- 
dard cosmology. In other words, we need to start with 
three spatial dimensions which at the time are much 
larger than the string scale (the entropy problem), and we 
can provide no explanation for why the current universe 
is so close to being spatially flat. Inflationary cosmology, 
in addition to providing a causal mechanism for structure 
formation, has the potential of solving the entropy and 
flatness problems. On the other hand, scalar field-driven 
inflationary cosmology is singular, whereas string gas cos- 
mology has the potential of providing a non-singular cos- 
mology. 
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